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Abstract 

The distribution of local charge excesses (DLC) in metallic alloys, previously obtained as a result 
of the analysis of order N electronic structure calculations, is derived from a variational principle. 
A phenomenological Charge Excess Functional (CEF) theory is obtained which is determined 
by three concentration dependent, material specific, parameters that can be obtained from ah 
initio calculations. The theory requires modest computational efforts and reproduces with an 
excellent accuracy the DLC and the electrostatic energies of ordered, substitutionally disordered 
or segregating metallic alloys and, hence, can be considered an efficient approach alternative to 
conventional electronic structure calculations. The substantial reduction of computing time opens 
new perspectives for the understanding of metallic systems and their mechanical properties. 
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The application of metallic alloys in a huge variety of high-tech areas, ranging from 
medical prosthesis [1] to jet engines Q|, requires a careful assessment of the mechanical 
properties of these materials. The determination of their phase diagrams is crucial to this 
end since the performances of alloys are heavily influenced by the various crystalline phases 
of which they are made 0, S, 0] . Recently, Zhao 7] introduced an innovative experimental 
method that allows for the rapid yet accurate assessment of alloy phase diagrams. However, 
theoretical methods with similar high-throughput periormaiices seem much beyond the status 
of art. 

Current theoretical approaches can be divided in two different classes. Theories in the 
first group are based on Landau's idea that ordering occurs in alloys due to some instability 
of the high-T solid solution phase (HTSSP) with respect to certain concentration fluctua- 
tions ^,0]. Such schemes often use the Coherent Potential Approximation (CPA) solution 
for the HTSSP as the reference state for a perturbation theory. They can be put in a 
very elegant form using a reciprocal space formalism, making contact with the concept of 
long-range order (LRO) and with the Fermi surfaces properties jlfll 11 1. In spite of many 



successes, these approaches suffered from the criticisms moved to the CPA which, in its 
standard implementation within the density functional theory (DFT), fails to account for 
the electrostatic energies in the HTSSP. Theories in the second class are based on Ising 
hamiltonians. The corresponding parameters can be extracted from the experiment, from a 
perturbative expansion of the CPA HTSSP solution or, as it is nowadays more common, 
from the total energies of various alloy configurations as obtained by DFT calculations j3| ■ 
In this case, accuracy is controlled by the cut-off length I assumed for Ising interactions. 
In metallic alloys LRO and short-range order (SRO) are entangled in such a way that the 
convergence with / is slow, while DFT calculations would be required for about A^! alloy 
configurations 14, Q|, with A^ ^ l^. In turn, DFT calculations for a_single configura- 



tion constitute a bottleneck: using the fastest available algorithms Il6l Il7|l the number of 

n 

floating-point operations required is proportional to A^, with huge prefactors Remark- 
able results have been obtained by a mixed approach that includes the summation of 
series both in the direct and the reciprocal spaces. However, in this approach no estimate 
of the size of the truncation errors is available. 

In the present Letter we shall show that the computational labour necessary to obtain 
total energies for any alloy conflguration can be greatly reduced. At variance of existing 
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simplified approaches [IJ], we shall maintain accuracies comparable to those of DFT 
calculations. 

The analysis 0, 0, of DFT calculations data in extended metallic systems indicates 
that the net charges at the crystal sites, g^, and the site Madelung potentials, Vi, are strongly 
correlated. More precisely |2Q,], within numerical errors, i) the DLC is continuous over a 
certain interval and ii) the pairs {qi,Vi), lie on straight lines, one for each alloying species. 
For future reference, the statement (ii) is conveniently rewritten as 

aiQi + Vi = ki (1) 

For a specified configuration of the binary alloy A^^Bc^, the coefficients in Eq. take 
the values a a and kA if the i-th site is occupied by an A atom or and ks otherwise. 
The results (i) and (ii), in the following referred to as the qV laws, have been numerically 



obtained within the Local Density Approximation 
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DA) and the muffin-tin or the atomic 
but not yet formally derived. Here 



sphere approximations for the crystal potential 
their validity shall be assumed as an "empirical" evidence. 

Accurate calculations of the alloy total energies must necessarily keep into account the 
above results. It has been shown recently that an isomorphous CPA model including local 
fields (CPA+LF) 



and a modified screened CPA approach j22| are able to catch the 
linear nature of the qV laws and compatible with (i). However, the same models are not 
able to self-consistently determine the DLC. Here we shall obtain the DLC from a variational 
principle in terms of one electronic degree of freedom for each atom, the local charge excess g,. 
In the resulting phenomenological Ginzburg-Landau theory, hereafter referred to as CEF, 
the site chemical occupations determine the qi at each site and the electrostatic energy. 
For all the (ordered or disordered) alloy configurations corresponding to the same mean 
concentration the theory is completely determined by three material specific parameters that 
can be obtained ab initio by supercell DFT calculations or by the CPA+LF theory. In the 



following the CEF theory shall be presented and 



;ested vs. order N Locally Self-consistent 



Multiple Scattering (LSMS) DFT calculations 2M fo'^ t)cc CuZn alloys. 

As in the case of LSMS calculations we study the binary alloy Ac^Bc^, ca+cb = 1, by 
the means of supercells of volume V containing N atoms with periodic boundary conditions. 
Each site can be occupied by an A or a B atom. The nuclear charge Zi and the volume [2^ 
J2i^i = V, are associated with each crystal site. Each alloy configuration is specified by 
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a set of occupation numbers [30|, , where = 1 if the i-th. site is occupied by a a atom 
or otherwise and I]jXf = Nca- The site charge excesses, qi = drp{f) — Zj, defined in 
terms of the electronic density p(r), satisfy the global electroneutrality constraint 

E?^ = o (2) 

i 

As it is verified within spherical approximations for the crystal potential and the LDA, we 
assume that the system total electronic energy consists of the sum of site-diagonal terms plus 

Em = J2ij Mijqiqj = ^J^iQiVi- The Madelung matrix elements 
23l . l2q , and the Madelung potentials are given by, 



a Madelung term j20|,|25| 
Mij are defined as usual 



V, = 2Y.M,,q, (3) 

j 

Unless otherwise stated we use atomic units in which = 2. 

We wish to develop a theory that determines the qi and incorporates the qV laws. For 
a particular alloy configuration, Eq. Q and the substitution of Eqs. Q in Eqs. ((T)) give 
X + 1 linear equations in the q^. This implies that not all the coefficients and ka can 
be independent. The CPA+LF model gives a hint for the missing relationship. In this 
theory j23| the quantities l/tta depend only on the mean alloy concentration and are viewed 
as the responses of impurity sites, embedded in the CPA 'mean' alloy, to a local field designed 
to simulate the Madelung potential. Furthermore, the CPA 'electroneutrality' condition for 
the zero-field charges, ba = k^/aa gives c^6a + 0^6^ = 0. These circumstances suggest that 
the theory we are elaborating should allow for a possible renormalization of the constants 
ka in different configurations corresponding to the same mean alloy concentration. 

To make further progresses, we consider the Ginzburg-Landau functional of the site charge 
excesses, 

ml^A = E |(^^ - ^^f + E Mim, -^^J2Q^ (4) 

i ij i 

where the Lagrange multiplier /i has been introduced to impose global electroneutrality. The 
minimization of Eq. (@]) with respect to the order parameter field {qi} and to fi gives the set 
of Euler-Lagrange equations constituted by Eq. (j2)) and by 

ai{qi - 6j) + 2 ^ Mijqj = (5) 
j 

Eqs. (0) are equivalent to Eqs. only when fi = 0. When /i 7^ 0, the renormalization 
H occurs to ensure global electroneutrality. Once the four constants Ua, ba have 



been obtained for a given alloy configuration, they also can be used for other configurations. 
As we shall see, the transferability of the parameters within fixed concentration ensembles 
is a peculiar strength of the present approach. 

Now, since Q has the dimension of an energy and contains the electrostatic contribution 
Em, we assume that, except but for an additive constant, its minimum value corresponds 
to the total electronic energy of the alloy configuration at hand. The quadratic terms in 
Eq. Q can be interpreted as energetic contributions due to local charge rearrangements and 
/i as the chemical potential ruling charge transfers. 

The explicit solution of the problem can be written in terms of A, the inverse of the 
matrix of elements 2Mij + J2a O'aX'^ 6ij. Straightforward calculations show that 

= [kA - ks) [(1 - y) E A.,^/ - y E (6) 

j j 

where Aap = J2ij Xi^AijX/ , and y = Y^a^Aa/ Yjap ^ap- The set of four constants found 
by LSMS or similar calculations for a specific configuration is then equivalent to the three 
constants in the CEF theory, a a-, clb and kA — ks, and to the chemical potential /i. The 
CEF formulation and solution are well defined for both ordered and disordered alloys. 

In the following we shall apply the CEF theory to bcc Cuo.50Zno.50 alloys and compare the 
results vs. LSMS calculations. We have selected several sample supercells described in Table 
I and designed to simulate random, partially ordered or segregated configurations. For each 
sample a set of the CEF parameters has been extracted by linear fits of the corresponding 
LSMS qV data: CEF calculations made using parameters from the n-th sample are indicated 
below as CEF-ra. 

A detailed comparison between LSMS and CEF-1 calculations for sample 1 is reported 
in Table II. The two sets of calculations present very small differences: 5 parts over 10^ for 
the mean values of the charges and of the Madelung potentials, 2 parts over 10^ for the 
Madelung energies, less than 1 per cent for the widths of the DLC. The absolute values of 
the differences Aqi = qf^^^^ — qf^^^^ are smaller than 0.005 electrons at any lattice site 
and not correlated with the chemical occupations. The mean square deviation between the 
two set of charges, ((Ag)^), is of the order of 10~^, i.e. it is comparable with the numerical 
errors in LSMS calculations. The resulting DLC's, plotted in Fig. (^, appear very similar. 
The main source of the tiny differences is that all the CEF charges hy construction satisfy 
the qV laws, while the same laws hold only approximately for LSMS calculations. 
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TABLE I: Characterization of the samples used and CEF parameters, a-cui ^Zn a-nd kcu-^Zn- All 
supercells correspond to bcc equiatomic CuZn alloys at the lattice constant a=5.5 a.u. SRO are 
the Warren-Cowley short-range order parameters for the first 3 neighbours shells. Samples 1-3 and 
samples 4-5 correspond, respectively, to random and partially ordered configurations, samples 6 
and 7 are, respectively, Cu/Zn and Cun.125Znn.s7R/Cun.s75Z1 in. 12 5 multilayers stacked along (001). 
LSMS calculations for sample 1 have been published in Ref. 



Samples 


N 




azn 


kCn-kzn 




SRO 




1 


1024 


1.84 


1.82 


0.29 


-0.004 


0.021 


0.010 


2 


256 


1.85 


1.83 


0.29 


0.078 


0.042 


0.021 


3 


256 


1.81 


1.80 


0.28 


0.031 


-0.005 


-0.021 


4 


256 


1.83 


1.83 


0.29 


0.184 


0.120 


-0.078 


5 


256 


1.90 


1.87 


0.30 


-0.309 


0.271 


0.141 


6 


256 


1.81 


1.83 


0.28 


0.750 


0.833 


0.667 


7 


256 


1.85 


1.85 


0.29 


0.414 


0.458 


0.375 



TABLE II: CEF-1 calculations for the bcc Cun.5nZnn.5n random alloy sample 1 compared with 
the LSMS results of Ref. I2JI. {q)cu and {V)cu are, respectively, the mean charges and Madelung 
potentials at the Cu sites, acu and azn the standard deviations of the DLC's for Cu and Zn sites, 
Em /N is the Madelung energy per atom and ((Ag)^) are the mean square deviations between CEF 
and LSMS charges. 

CEF-1 LSMS 
{q)cu 0.099787 0.099783 

{V)cu -0.038197 -0.038188 

acu 0.02507 0.02523 

azn 0.02801 0.02814 

Em/N (mRy) -2.552 -2.557 

((Ag)2) 2.7 10-6 
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FIG. 1: Cu (light histogram) and Zn (dark histogram) calculated DLC for the bcc Cuo.soZno, 



random alloy sample 1. Top frame: LSMS (Ref. 
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2J); lower frame: CEF-1. 



Specific tests have been designed to check the transferability of the CEF parameters 
extracted from one sample to the others. The mean square deviations of the charges obtained 
from CEF-n and LSMS ((Ag)^) (Table III) clearly show that there is no appreciable loss 
of accuracy when the CEF parameters extracted from random samples are used in ordered, 



partially ordered or segregated samples and vice versa. Calculations j27| for CuPd alloys 
support the same conclusion. The transferability of the CEF parameters is a very remarkable 
result and implies that the theory is generally applicable to metallic alloys, no matters 
whether they are ordered, disordered or segregated. Furthermore, the very high accuracy 
obtained for the Madelung energies (Table IV) indicates that the theory can describe very 
carefully the electrostatic contributions to the energetics of ordering phenomena. 

We conclude this Letter with remarks about interesting aspects of the CEF model in view 
of possible future applications. 

The CEF operates a coarse graining over the electronic degrees of freedom that are re- 
duced to one for each atom, the local excess of charge. This notwithstanding, the theory 
carefully reproduces the DLC and the energetics of metallic alloys in any ordering status. 
This has been possible because the Madelung potentials, through Eq. 0, weight as appropri- 
ate the long ranged effects of the occupations of all the crystal sites. Such a renormalization 
of the interactions holds exactly within CPA-based theories like the CPA+LF or the 
PCPA |2^, where any site diagonal property is a unique function of the Madelung potential 
Vi and the nuclear charge Zi at the same site. This uniqueness does not hold j22| for more 
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TABLE III: Charge mean square deviations ((Ag)^) x 10^ between CEF-n and LSMS calculations. 



Columns identify different samples. 

1 2 3 4 5 6 7 

CEF-1 3 3 3 2 7 0.4 3 

CEF-2 3 3 3 2 7 0.4 3 

CEF-3 3 2 3 1 9 0.5 3 

CEF-4 3 3 3 2 7 0.4 2 

CEF-5 4 4 5 3 4 0.5 3 

CEF-6 2 3 3 2 6 0.5 2 

CEF-7 4 4 4 2 5 0.4 3 



TABLE IV: Madelung energies per atom (in mRyd units) from LSMS and CEF calculations. 
Columns identify different samples. 





1 


2 


3 


4 


5 


6 


7 


LSMS 


-2.56 


-2.15 


-2.32 


-1.42 


-3.70 


0.20 


-1.02 


CEF-1 


-2.55 


-2.11 


-2.30 


-1.40 


-3.82 


0.21 


-1.03 


CEF-2 


-2.54 


-2.10 


-2.29 


-1.40 


-3.80 


0.22 


-1.02 


CEF-3 


-2.57 


-2.12 


-2.31 


-1.41 


-3.85 


0.20 


-1.05 


CEF-4 


-2.54 


-2.10 


-2.29 


-1.40 


-3.81 


0.21 


-1.03 


CEF-5 


-2.49 


-2.05 


-2.24 


-1.37 


-3.71 


0.25 


-0.98 


CEF-6 


-2.52 


-2.08 


-2.27 


-1.39 


-3.78 


0.20 


-1.02 


CEF-7 


-2.50 


-2.06 


-2.25 


-1.37 


-3.73 


0.23 


-1.00 



exact approaches, where some residual dependence on the site nearest neighbours environ- 
ment is expected for. Nevertheless, the fact that CPA-based theories accurately accounts 



for the spectral properties of metallic alloys 



21[ and the quantitative agreement with LSMS 



231 and |28|, suggest that the errors in- 



calculations, in the present work as well as in Refs. 
troduced by neglecting the nearest neighbours effects not already conveyed by the Vi are 
comparable with numerical errors in DFT calculations. We mention that, in a very different 



context, precedents of this idea of a coarse graining over quantum degrees of freedom can 
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be found in the concepts of chemical valence and of electronegativity. 

The calculations presented in this Letter require floating-point operations in order to 
obtain A by conventional linear algebra algorithms. Thus, for = 1000, the CEF is about 
10^ times faster than LSMS We are convinced that, in the next future, such a large 
computational speed up and the mentioned transferability of the CEF parameters will make 
possible the development of accurate ab initio techniques for the investigation of ordering 
phenomena and the calculation of phase diagrams in metallic alloys. 

We thank J.S. Faulkner who made the data of Ref. 
acknowledge discussions with E.S. Giuliano. 
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